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ABSTRACT 


Vertical motion resulting from forced ascent or descent over 
sloping terrain is the primary concern of this investigation. In 
addition, a method of computing precipitation is adapted to estimate 
the contribution made to the total precipitation by the terrain-induced 
vertical velocities. In mountainous regions it is assumed that hor- 
izontal flow at levels above the surface is non-divergent and may be 
described by the stream-function, A direct method of stream-function 
calculation is used to yield the stream-function at the locations where 
the wind is measured, An objective analysis, which utilizes a con- 
tinuous distance weighting scheme, extends these values to a grid 
oriented such that its ordinate parallels the Coastal Mountains and the 
Alberta Range of the Rocky Mountains in southern British Columbia and 
Alberta. 

A scalar objective analysis is proposed to compute the hori- 
zontal velocity components relative to the grid from the more abundant 
reports at the earth's surface, This analysis is also used to estimate 
values of temperature and dewpoint temperature at the grid points from 
station reports, In the latter case, the analysis was applied not only 
to the surface, but also to the 850-mb, 700-mb and 500-mb levels where 
station density allows only a first approximation, 

The terrain-induced vertical velocity is calculated from the 
equation: 


Ww = Vb . Vp, H 9 
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where Vy, is the horizontal wind and VY )H, the slope of the terrain is 
decreased linearly with height to zero at 3 km above the terrain, 

The orographic precipitation is computed using the vertical 
velocity, temperature and dewpoint temperature fields at a number of 


levels and compared to surface measurements of precipitation, 
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CHAPTER 1 
INTRODUCTION ~ 


All extensive mountain systems around the world influence 
weather from the planetary scale down to the microscale, The inter- 
actions between topography and the atmosphere are exceedingly complex; 
attempts to study these relationships have met with varying degrees of 
success, but as yet a complete understanding is lacking. Part of the 
problem arises from the overlap of the various meteorological scales, 

In general the terrain may be highly irregular, resulting in rapid 
changes of forces exerted on a particular parcel of air. Also to be 
considered are the overall orientation of the flow with respect to the 
obstruction, the stability of the air mass, vertical and horizontal 
Wind shear, and moisture content to name only a few, 

An evaluation by Newton (1971) of the torques exerted by large 
mountain ranges on the atmosphere suggests that this factor may be com- 
parable to surface friction in the balance of global angular momentum, 
On the synoptic scale, Pettersen (1956a) in a study of cyclogenesis over 
the Northern Hemisphere noted that the lee of extensive mountain 
barriers was a major ache for the formation of lows, A more recent 
study of this same phenomenon by Chung (1972) considered the area to the 
lee of the Canadian Rocky Mountains. By reducing the size of his samp- 
ling squares, Chung found two additional major centres of cyclogenesis 


that were not revealed by Pettersen's investigation. The last example 
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is typical of the problem confronting meteorologists. A further re- 
duction of the mesh size would in all probability reveal an even more 
involved pattern, if data were available, 

McClain (1960) found that descending motion to the lee of the 
Western Cordillera of North America was not in itself sufficient to 
generate major cyclonic storms but, when coupled with other factors, 
could result in rapid and dramatic development. Katabatic winds, moun- 
tain waves and other mesoscale phenomena usually attributed to the 
proximity of marked relief are also of great importance but influence a 
much smaller area, 

One of the most generally known aspects of mountain meteor- 
ology is the enhancement of precipitation along the windward slope of 
any orographic obstacle. Sawyer (1956) gives a brief but comprehensive 
summary of the problems involved in estimating orographic precipitation 
from physical and dynamical principles. Reinelt (1970) points out that 
while in some cases the rain shadow effect may predominate, a change in 
the circulation pattern will make the same area subject to terrain- 
induced rain. 

The present study was based on the idea that air encountering 
a large mountain range will be forced to ascend the windward slope and 
descend on the lee. A model proposed by Charette (1971) to calculate 
the terrain~induced vertical velocities at the 850-millibar level was 
modified so that less input data were required, This allowed his 
approach to be extended to higher levels where the density of data dec- 
reased rapidly. 

Chapter 2 contains a description of the grid and the vertical 


velocity calculation, The method of determining the horizontal wind 
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field at discrete grid points is outlined in Chapter 3. A quantitative 


precipitation scheme for evaluating the vertical velocities computed by 


the model is given in Chapter 4. The results of the calculations are 


presented in Chapter 5. 
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CHAPTER 2 
OROGRAPHICALLY-INDUCED VERTICAL VELOCITIES 


21 Introduction 

McClain (1960) observed that the Western Cordillera of North 
America constitutes an essentially continuous barrier from 25°N to 
60°N for air below 1300 to 1500 meters. The 850-mb level is usually 
found near the upper limit of this barrier and must therefore reflect 
its influence. When presented with so great a north-south extent, most 
of the air below 1300 meters must surmount the obstacle and-descend on 
the lee slope, The method used to calculate the resultant vertical 
motions was based on Charette'’s (1971) study, Some modifications were 
made to the general procedure and these are described in the following 


sections. 


Zen ine Grid 

Charette extracted terrain heights at discrete points from 
topographic maps of scale 1:250,000. These maps were part of the top- 
ographic series produced and printed in Canada by the Surveys and 
Mapping Branch, Department of Energy, Mines and Resources, Government 
of Canada and in the United States by the Interior Geological Survey, 
Washington, D.C... The location of these points was a function of a 
polar stereographic projection true at 60°N, A brief description of the 
equations relating a point on the earth's surface to its projection onto 


some image plane is given by Saucier (1955) and in more detail for the 
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particular case of a polar stereographic secant projection by Charette 


(1971). 
The relation may be determined from: 
r = a(1 + eos (85) (tan (2) (2-1) 
2 (27 coe (eon) mle sin (5) 
Oo = (22) 
(is "cos" (6) (1 + sin (0)) 
where 
r = radius of a latitude circle on the polar 
stereographic projection 
a = mean radius of the earth 
> = standard latitude, in this case 60°N 
@ = latitude of point in question 


96 = 1/2 Se 0 
e = 1/2-% 
o = image scale which gives the distortion of 
distances introduced at the latitude 0 by 
the representation of the earth's surface on 
a polar stereographic projection with stan- 
dard latitude Qo. 
Only equation (2-1) is required to compute the location of the grid 
points, The grid most commonly used by numerical weather prediction 
models has a grid length of 381 kilometers at 60°N on a polar stereo- 
graphic projection, This grid size remains constant on the projection 
but the actual distances between grid points varies according to the image 
scale oe In this investigation a grid length equal] to one-eighth of 
381 kilometers was chosen as this was the same used by Charette in the 


extraction of terrain heights, The ordinate of the grid was aligned 


with longitude 80°W. 
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FIGURE (2-1) Calculation of the grid point coordinates from a polar 
stereographic secant projection at 60°N, 


As illustrated in Figure (2-1) the distance r from the pole to 
a grid point on the projection plane may be computed in terms of 


one-eighth grid length? intervals. 


r= (xa)? + (ya) 1/2 


where d is the grid length or 


72 
r = a[x2 + y2] / (2-3) 
Note that X and Y represent the number of grid lengths along the absci- 


ssa and ordinate, respectively. Values of increasing Y are measured 


lfhe one-eighth grid length will henceforth be referred to as 
the grid length, unless otherwise stated, 
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downward, and increasing X to the left. By substituting this expression 


for r into equation (2-1) it is possible to solve for the latitude on 


the earth's surface: 


2 We 
ie ete reese [ellen (2-4) 
2 al leticos+ (0. 


With \ defined as in Figure (2-1) the longitude of a gridpoint may be 


determined from equation (2-5): 


»~ 
u 


80° + A~ 


fl 


80° + arctan (X/Y) 


The portion of the map used in the analysis encompassed all of 


British Columbia, all of Alberta with the exception of the northeast 


corner, and parts of the adjacent regions of Montana, Idaho, Washington, 


Oregon, Northwest Territories, Yukon, Saskatchewan and Alaska (see 


Figure (2-2)). 
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FIGURE (2-2) Grid used for orographic vertical velocity model, 
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Ze3, Ihe Height Profile 

Once the coordinates in latitude and longitude of each grid 
point were calculated the tedious but necessary extraction of heights 
from topographic charts was performed by Charette. He relates that, on 
the 1:250,000 scale, contours ranged from every 100 feet over relatively 
level regions to 500 feet in areas of high relief, Thus, the accuracy 
of the terrain heights varied as the ruggedness of the topography. 

After the heights were procured he performed a smoothing 
operation on the data at every second grid point. The operator used is 


given by equation (2-6) with reference to Figure (2-3): 
H(I,J) = [e(nz tho +hg +hy +hs + hg + ho + hg)]* 0.6 


a NGS) bs On (2-6) 


O1/4 grid points 
e 1/8 grid points 


1/8 grid length 
i | ant 


FIGURE (2-3) System of grid points used in smoothing operation. The 
height profile is retained at every second grid point following applica- 
tion of smoothing operator, 

The effect of the operator combined with retention of the 
smoothed profile at alternate grid points provides a simplified but dis- 


torted height field. A portion of the original data field analyzed 
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before smoothing is shown in Figure (2-4). For comparison, the same 
area using the smoothed values is shown in Figure (2-5). 

The reason for smoothing the field is to eliminate irregul- 
arities which, it is argued, have no bearing on the synoptic scale; how- 
ever, the reduction in height of major ridges, and the associated relax- 
ation of height gradients can be severe and excessive, In this study, 
it was decided, therefore, to retain the original terrain heights with 
the realization that a number of "“unrepresentative” values might be 
present, This problem exists whenever any attempt is made to represent 
a continuous field by discrete, equally-spaced values, 

The justification for not smoothing may be found in the mean 
value theorem from calculus. In this investigation, the important 
quantity was taken to be the slope of the underlying terrain. From the 
theorem it can be stated that at least one point between any two grid 
points must have a slope equal to the height gradient between the grid 
points, wrovided a first,/Cerivative exists\at.all the intervening 
positions, Therefore, the data were utilized without subsequent smoo- 
thing. The only operation performed on the data was a conversion from 
feet to meters, An analysis of the complete area used in the model is 
given in Figure (2-6), With reference to the grid shown in Figure (2-2), 
the analysis extends from rows 28 to 85 and columns 40 to 68 inclusive. 
Some distortion of shape and lack of prominent peaks results from the 
discrete representation, but it was felt this approach resembled reality 


more closely than a smoothed profile. 
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FIGURE (2-11) Analysis of terrain heights extracted by Charette, 
heights are labelled in thousands of feet above sea level. 
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FIGURE (2-5) Analysis of smoothed terrain showing the same area as in 
Figure (2-4). The heights are labelled in thousands of feet above sea 
level. 
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FIGURE (2-6) Height profile of that portion of the Western Cordillera 
used in this model. The heights are expressed in hundreds of meters 


above sea level. 
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2.4 The Vertical Velocity Calculation 

One of the most difficult quantities to evaluate in meteor- 
ology is vertical velocity. The problem stems from its small magnitude 
so that it is difficult to measure directly and accurately in the free 
atmosphere. Most studies of ascending or descending atmospheric motions 
determine the vertical velocity indirectly from other more easily 
measured parameters and theoretical principles, 

Three scales of motion are commonly associated with the form- 
ation of precipitation, as mentioned by Reinelt (1970). Vertical vel- 
ocities of a few centimeters per second that are maintained for at 
least several hours are generally attributed to synoptic systems such as 
extratropical cyclones ( Mason, 1962, p.9 ). Near fronts this figure 
may increase to 10 cm/sec, but this is considered the upper Limit -Or 
organized vertical netion over a large area, At the other end of the 
scale are convective motions which may exceed 20 meters per second in 
large thunderstorms, but are more usually of the order of 5 m/sec. An 
individual convective cell seldom exceed a few kilometers in horizontal 
extent or half an hour in duration. The third major range of vertical 
velocities intermediary to the others is associated with forced ascent 
or descent due to the presence of orographic obstacles. 

The last case is the one studied in this investigation. Ver- 
tical velocities greater than 40 em/sec on this scale are considered 
rare as this requires strong horizontal winds impinging on an extremely 
steep slope. Such slopes are usually of small geographical extent, so 
that it is probable that air is diverted to adjacent regions, Along a 
mountain ridge, vertical motion of 5 to 15 cm/sec is more in order, 


Over gently sloping terrain prolonged ascent at less than 1 em/sec for 
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a fetch of hundreds or even thousands of kilometers has also been 
observed, Of these three situations, the second is thought to be of 


prime importance in the mountainous regions of the world. 


FIGURE (2-7) Cross-section of a ridge of land illustrating orographic 
ascent, 


15 


A parcel forced to ascend a mountain (Figure (2-7)) is assumed 


to follow the shape of the underlying terrain, By a method described 
Chapter 3, the horizontal velocity Vp, is obtained at the centre of eac 
grid square at a number of levels. At the surface the vertical veloci 


was computed from equation (2-7): 


W = Vn a VpH (2-7) 


where Ww = vertical velocity component and 


Vy = horizontal wind 
v4 = yao 4 oo the slope of the terrain. V)H 
ox oy 


was evaluated at the centre of the grid squares, The height at point 


A (see Figure (2-8)) was the average computed from: 
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Ha = 3(H(I,J) +H(I +1, J)) (2-8) 


and, similarly, the heights at points B, C and D. In actuality this 


step was incorporated in the expressions for the gradients given below. 


oH : 
ee) se eC 11) Bead tad) 
(2-9) 
oH 0 
Same = ayy RB em Gay le) Sat cra ape (Ger hears ay) 
(2-10) 
where go = image scale from equation (22) 
d = grid length in meters 
In component form, equation (2-7) may be written: 
, 3H aH 
= y—— + vy— no bel 
W vere ot ag (2-11) 
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FIGURE (2-8) Evaluation of slope of terrain at point E, the centre of 
a grid square, 


The influence of the orography was assumed to decrease with 


increasing height. McClain (1960) stated that no evidence existed at 
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that time to indicate that the synoptically important orographic motion 
extended more than 3 kilometers above the terrain, While the scale 
considered in this study was much smaller than synoptic, this was used 
as a criterion for the upper limit of orographic vertical motion, Mountain 
waves are ignored, | 

The exact nature of the decrease has not yet been determined 
for a general case although various proposals have been made by 
Smagorinsky and Collins (1955), Estoque (1957) and numerous others. 
Although there was no fundamental physical reason for doing so, it was 
decided to follow Charette's lead in this matter and allow the effective 
slope of the topography to decrease linearly with height to zero at 
3 kilometers above the terrain, Let 8H represent the distance in meters 
between the level under consideration and the terrain at Peeeeentre ofa 
grid square, Then the vertical velocity at that level due to the oro- 


graphy may be expressed as: 
eee) yp (2-12) 


In practice, the correction is applied following the calculation of w as 
given in equation (2-7). The sole justification for this approach is 
its simplicity, although later evaluation indicated that it was useful 
as a first approximation. More elaborate methods await further invest- 
igation. 

The major difference between Charette’s calculations and those 
outlined above was the reduction of the grid length used to calculate 
the slope of the topography. By centering the calculations at the middle 


of the grid squares it was possible to eliminate consideration of the 
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wind direction when computing the gradients. Downward motion may now 


be overestimated in the case of narrow valleys where air tends to stag- 


nate but this was not thought to be important in the summer and autumn 


situations studied, 


CHAPTER 3 


STREAM-FUNCTION ANALYSIS AND SCALAR OBJECTIVE ANALYSIS 


Dal Lneroduc Tien 

In order to proceed with the vertical velocity calculations 
mentioned in the previous chapter it is necessary to compute the hori- 
zontal wind components at the centre of the grid squares. Unless reports 
are numerous, an interpolation scheme between station reports results 
in a poor estimate. A much more powerful approach was used by adapting 
a direct method of stream-function calculations from actual wind reports 
as proposed by Endlich and Mancuso (1964). This iterative procedure 
permits evaluation of the stream function at the stations, In a later 
report, Endlich and Mancuso (1965) included an objective analysis for 
the computation of the stream-function at grid points from those deter- 
mined at the stations, 

Charette (1971) applied the method of Endlich and Mancuso with 
only minor modifications to an area slightly larger than that indicated 
in Figure (2-2). However, some major changes were introduced into this 
scheme in the present study. 

ithe Erebe wine field was derived from a scalar objective 
analysis using a correction method which bears some similarity to that 
proposed by Cressman (1959). This essentially amounts to an elaborate 
interpolation scheme but, at the earth's surface, it was felt that reports 


were sufficiently abundant to allow this. It was found that with a 


7 = 7 i - 
—_ : i =e A Tt Sarre anth ise 
‘ ‘ > @ a] uid TEI Ae iy ‘iva c TARA heey at? yas? 
7 


: a 
‘ 
avs = 
& pas. 
ors 
a 
; : : aad « o A Z 
; re. = ; ~@ ep Ty oa i 
im@se.ht 4 ome! is f f is Ms 
iMivaw a . 
a 
[ 
p 
* el yaeriea ' > 
. it migitAs ‘ote Satna 
- ‘ ! "1 «4 j o 7 4 1a ceo iso. < wad 7 
ia 
4 y —_ ’ 
os 
’ 1 - i io e l 5 
G 
@ 
a3 z 4 4 j 7 Ps 7 
cal 
. 
qv , 
j a * 
Ris 
: = 
i a 
- ‘ ae a . 
I 
= 7” 
i a , oa } po a jm Blea 


refine 
- fh . 
= 7 a> . ' 
#7 ¥ 7 a yo pa \ q vo) Pt eee 
ms 7 — 


: = - 
eid? . arn? Pe oh aes : PPh 3.4 (i Serart Ps fat 
nt | . iba 
+ ; ‘ : i 7 


4 farts 
a 


tk 


20 


minimum of modification this same procedure could be applied to the 
temperature and dewpoint fields at ground level, The scalar objective 
analysis is also used at upper levels to estimate the temperature and 
dewpoint at grid points with the realization that the derived fields 


may be in error in data sparse regions, 


3.2 The Stream-Function: A Brief Review 


Consider flow parallel to the x-y plane for an incompressible 
fluid. It can then be shown (see e.g. Whitaker (1968)) that the mass 
flow rate across an arbitrary curve between two points is independent of 
the path chosen, That is Ven ds is an exact differential of some func- 
tion , where v is the velocity vector along the curve, n Loe a) Unit 
Telit normal to the curve, and ds, an elemental distance along the curve, 


This function is known as the stream-function, and the differential dy 


expressed as a function of x-y coordinates is: 


dy 


(SBoax + (SB)ay (3-1) 


but 


dy ¥Vnds = un, ds +v ny ds (3-2) 


Aa 


> > 
where u, v and ny, ny are the components of V and n. Let t be a unit 
vector tangent to the curve joining the points mentioned above, It is 


easily shown that: 


— eee = ax + dy = 
Ont xe tay (a3) 
and therefore 
Nx aa and ny ae (3-4) 
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The signs in expression (3-4) could be reversed and still satisfy 
equation (3-3), but the signs in (3-4) are retained in the subsequent 
steps. By equating equations (3-1) and (3-2) and substituting for nx 


and ny from (3-4), 


ow ow 
te Oe mavens Os (3-5) 
As equation (3-5) applies for any curve with the implication that dx 
and dy are not identically zero in every case, it follows that: 


ow 
Ox 


uio= and v = ‘ (3-6) 


yall 
— 
Hess (1959) defines a streamline in the following manner: 
A streamline is a line to which, at any given 
instant, are tangent all the velocity vectors of the 
points through which the line passes; thus the flow 
is along the streamlines at any given moment, 


Along a streamline the velocity vector is \¥| t and the difference in 


the stream-function at two points A and B: 


3 eee 
vgp-¥, = Sav = J Ivitnas = 0 (3-7) 
Thus, at a given instant, a line joining points of a constant stream- 
function is a streamline. 

In the atmosphere, flow may be treated as incompressible if 
the divergence is zero or vanishingly small. Like vertical velocity, 
divergence in the atmosphere cannot be measured directly because of its 
small magnitude and the errors in wind measurements from pilot balloon 
or radiosonde ascents, The next section is based on the assumption that 


the divergence is zero so that the two dimensional flow is completely 
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described by the stream-function which accounts for vorticity, trans- 
lation and deformation, Application of this assumption to mountainous 
regions may result in a grossly simplified picture of the true flow, 

but it is believed that such a treatment is considerably better than the 


assumption of a geostrophic balance, 


3.3 Stream-Function Calculation at Stations 
Substitution of equation (3-6) into (3-1) yields the stream- 


function differential dy between two nearby points: 


dv = vdx -udy (5-8) 
Orn, 10 tinive difference forn, 
at Eee GN, - (3-9) 
The stream-function has dimensions (L“/T). It is convenient to define a 


modified stream function: 


* fa 
= eee -10 
Se tr) (3-10) 
where f = Coriolis parameter, The modified stream-function now has 


the same dimensions (L*/T2) as geopotential height, This manipulation 
allows comparison with the original height field which is also expressed 


in geopotential units, Therefore equation (3-9) becomes: 


ieee ees, = eB 
AY . = 5.5 (vAx = udAy) (3-10a) 


The difference in the modified stream-function between a station o and 
a nearby station i is: 
x k * 
AVs ae Vo = ai e 
Thus, the modified stream-function at the station o may be estimated from: 


* * * 
re weet tc 


- 7 2 : = - 
don bp y / eae nt Om: = (ae — 
: we >i ; ; = van 


yroa 


: a a 
i e ae NG | = - 


22 


* 
In their earlier paper, Endlich and Mancuso (1964) evaluated Avs; from 


the expression: 


ey oo reas (Yo f Vike a ne ue ye ve) 
i 9.8 iy Sree ewes 


Cebailaby 
where g has been replaced by 9,8 in agreement with equation (3-10). 
In the free atmosphere, flow is very closely geostrophic and is given 


by the vector product: 
Y= 22kx wz, (3-12) 


where Z is the height of a constant pressure surface in geopotential 
meters (gpm) and k is a unit vector normal to the horizontal plane. The 


stream-function also determines the wind from: 


“ n x 
V=kxW = os kx W (ais) 


In middle latitudes, the geostrophic balance is a good approximation to 


the actual wind field so thats 


In their 1965 version of the stream-function computation, Endlich and 
Mancuso used a weighted combination of winds and heights given in 


equation (3-15): 


een os Pty | v, +v3)(x, - x;) 
Avs Creve WT [( fe) i) (Xp i 
- (ug t+ uy)(¥9 - yi)] + .W3(Zy - 24) (3-15) 


where w} and W3 are weights with empirically determined values of 0.45 


and 0,10 respectively. 
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Extensive testing of both equations (3-11) and (3-15) revealed 
rather large discrepencies between the input height field and the final 
field of the modified stream-function. Changes as large as 60 gpm were 
observed on occasion, and, more frequently adjustments greater than 
30 gpm. Convergence within a prescribed limit of 0.5 gpm was much faster 
using equation (3-15), but some major differences appeared between the 
two modified stream-function fields when the same data were used for 
input, 

Despite the disadvantage of requiring more computer time, 
equation (3-11) was used in the final model, as the assumption that the 
flow was geostrophic in mountainous regions, needed to justify equation 
(3-15), was apparently violated. Charette used equation (3-15) in his 
model. 

More than one estimate of Vo is possible when the number of 
neighbouring stations is high. If there are n such estimates, then an 


averaged estimate might bes 


n 

co “ 

Vo = + 2 Coy + avy) (3-16) 

n 
Las 1 
A procedure suggested by Endlich and Mancuso (and followed by Charette) 
* 
was to compute two such averages, One average (Vo) 1 is obtained from 
all the stations within 300 miles of the station o, and a second aver- 
* 

aged estimate (Vo)o from stations located between 300 and 550 miles from 


station o, These two estimates are combined as: 


* * * 
Yo = Ar(¥o)y + Aalod2 » (3-17) 
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where Ay and A> are weighting factors assigned values 0.8 and 0.2, 
justified by the statement that: 
the observations are weighted more or less 
inversely with distance from the point of 
interest se. . 
(Endlich and Mancuso, 1964) 
A search for a better weighting scheme was undertaken, since 


the one given above is not only awkward to program, but is also discon- 


tinuous and arbitrary. An alternative which was tested is simply: 


(Ry BEAR 
A= 1-(r/p)* = 2S forr<R 
R 
(3-18) 
A = 0 agora Gal oe FR 
where R = a set radius of influence 
r = separation between two stations. 


It was found that although a generally smooth pattern was attained, in 
those areas where there was greater station density, the calculated 
stream-function was not particularly sensitive to an individual report, 
The wind field from such an analysis did not agree well with a part- 
icular wind report but was an average of the winds at surrounding stat- 


ions. Cressman used: 
es 2 2 
ea EY = Roos for rie R =] 
EaOAIE poriee (3-19) 


which gives slightly less weight to the stations in the range O<r<R 
(see Figure 3-1) than equation (3-18), A linear decrease with distance 


is given by: 


A = 1 (r/R) = oe for r <R (3-20) 
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WEIGHING FACTOR A; VERSUS STATION SEPARATION f; 


RADIUS OF INFLUENCE R=550 km 


I- (t/7R)* 


MANCUSO WEIGHTING SCHEME 


08 


0,7 


I= (r/RIC 
i+(r /R) 
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FIGURE (3-1) A comparison of various weighting functions dependent upon 
distance of separation. 


ar 


Neither of the last two weighting functions was considered significantly 
different from expression (3-18). A weighting factor which is close 
to unity when the station separation is small and decreases to near 


zero much more rapidly than any of above with increasing station separ- 


ation is: 


Amiaumexst (33/o)6 (3-21) 
where c is a constant which may be set to any value greater than zero, 
In this study the value of c was set equal to 200 and r was expressed in 
kilometers. 

Equation (3-21) determines the relative weights of the estim- 
ates of the stream-function at a particular station from its neighbours 
with a radius of influence of 600 kilometers, A weighted mean is given 


by: 


AG( V3 + Av; ) 
age dpm Nt Na eae (3-22) 


iT} feat fet 


Reports of wind speed and direction, and geopotential height 
were extracted from plotted and analyzed synoptic charts at &50 mb, 
700 mb and 500 mb. Pilot balloon reports did not have the heights of 
the pressure surface at their locations so that bogus data were inserted 
for the heights. In the first scan, these abnormal values were recog- 
nized and not used as initial guesses of the modified stream-function, 
Heights from radiosonde reports were used to establish the estimate of 
the modified stream-function at a pilot balloon station by the stream- 


function calculation, The new estimate of Vo was immediately substi- 
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tuted for the former estimate in much the same manner as a sequential 
relaxation. This had two main advantages over updating the field simul- 


taneously after a scan of all stations was completed: 


(a) convergence of v at all stations to within .5 gpm 
was much faster. 


(b) computer storage requirements were reduced. 


On any particular scan a record is kept of the maximum 
alteration between the present and previous estimates of ‘ . Lretnis 
maximum value is greater than .5 gpm another complete pass of all 
stations is made, As soon as the maximum alteration drops below this 
limit, the computation is terminated. In no instance did the field fail 
to converge, although it was noted that the magnitude of the alterations 
decreased very slowly from 1 gpm to .5 gpm, particularly with a large 
number of stations. Convergence to within 1 gpm usually took from 6 to 
8 iterations independent of the number of stations, It is possible 
that machine truncation error was responsible for the effects mentioned 
above, From experience, it is now Pelt that 1 gpm would be a more 
prectical limit. 

Another feature of interest was that a slightly faster conv- 
ergence resulted when the stations were ordered in such a way that the 
computations began with stations along the outer edge of the map and 
"spiralled" in toward the centre in the initial testing, Data were there- 
fore ordered in this fashion in all the cases investigated in Chapter 5. 

No adjustment was made to the final field as was done by 


Charette or Endlich and Mancuso, because the station density was not 


lm ge DS BR bsec f° 


a 


uniform in all regions. Charette's adjustment method is open to 
. ’ 1 > * 
objection on the grounds that it also changes the gradient of the over 


the entire map. 


3.4 Objective Analysis of Modified Stream-Function 


x 
Once the values of ~» are known at the stations it is necessary 
x 
to estimate » at the nearby grid points, Let 2 be a grid point near 


ra 
station i, Then yy at & may be calculated to be: 


* * ee 

Vo = Vt Ab (523) 
where 

* fpr Tae iva pee Us yy my) 

Av, = a sca nn en (3-24) 


If there are n such stations then: 


n 


i_™m 


x x 
Ai( v3 + Adz) 
1 


las (3-25) 
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where A; is a distance-dependent weighting factor as given in equation 
(3-21), Again this differs from the weighting scheme used by Endlich 
and Mancuso (1965), but it is considerably easier to implement and not 
discontinuous, 

At the sacrifice of using increased storage, the calculations 
are performed from the point of view of a given observing station. From 
the station's latitude and longitude it is possible to determine the 
station location in grid coordinates. A radius of influence of 600 km 
is used. Many of the stations used in the calculation lie outside 


the grid and may only influence a few of the outer rows or columns, The 
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range of indices of the rows and columns which include all the grid 
points within 600 kilometers of the station can be computed, Points near 
the corners of this region are not within the 600 km criterion, but as r 
is calculated for the weighting function anyway, these points are easily 
excluded, A weighted estimate and a weight are computed for each grid 
point within 600 km. Then the next station is considered. Each new 
weighted estimate is added to the sum of previous estimates at the grid 
point, as are the weights. This requires two matrices equal to the size 
of the grid. Once the stations are exhausted, the weighted mean is 
calculated at each grid point from equation (3-25). 

The modified stream-function may be used to evaluate the u and 


v components of the wind with respect to the grid from: 


- - 28 2 (3-26a) 
uy = v2 ay 3 a 
x 


If one replaces H by if in Figure (2-8) the following finite difference 
approximations may be used to evaluate the wind components at the centre 


of each grid square: 


x x * 
w= - 285 (H(I) + HET $1) -9 (Ets) - VD +4, 0 +1) 
v uel 
(3-27) 
and 
. * 
vy 2:82 ($(1,5) + ¥(r +1) - 9, 941) - $41,042) 
(3-28) 
where o = image scale 
d = grid length at 60°N 
f = coriolis parameter 


) is also-used as an eStimate of the height of the constant pressure 
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surface. If the mean value of  W at the centre of the grid square as 
; * 

determined by an average of w at the corners is less than H, aUDen 11 

and v are set equal to zero since the pressure surface would be below 


the terrain and no flow exists. 


Go) sl neiocalar Objective Analysis 


A means was sought to estimate the wind field at the surface, 
The assumptions of a geostrophic balance are not valid because of the 
proximity of the surface which exerts a frictional force on an air 
parcel, Divergence at this level is large so that the stream-function 
cannot be used to represent the total wind field. As a rough approx- 
imation the wind at the stations was expressed in component form with 
respect to the grid. IfS is the wind speed and D is the figeion of 
the wind in degrees measured clockwise from true north, then the comp=- 


onents ares 


u S cos(350 -D = A) 


(3-29) 


Vv S sin(350 -D =- jd) 
where } = longitude of the station.- 

On the first scan, a radius of influence of 700 km is used, 
The method of applying the estimates of u and v at the grid points is 
identical to that given in the previous section =e the exception that 


* 
equation (3-30) lacks a term analogous to AY, Let X be any variable 


for which a value is known at the station. Then: 


Xx, = (3-30) 


Isee Appendix A for details of derivation. 
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where A; = weighting factor as before, As may be expected, the first 
pass yields a very smoothed representation, Therefore, three more 
scans are made reducing the radius of influence by 200 km each time, On 


each scan after the first a correction term Cc) is computed from: 


2 ° 
ft . Ay( X; = Xq) 
al = 
Creo n (3-31) 
Bo Ay 
i=1 
and ronesin, CNL (3-32) 
Q Xo Q oe) 


If only one station influences a grid point then: 
Co = ASC Xe =X) (3-33) 


The final scan has a radius of influence of 100 km which is sufficient 
to give reasonable resolution. It might be practical to change the value 
of c in equation (3-21) but such a refinement was not investigated, 

Originally this method was designed to derive an initial guess 
field for a scalar objective analysis developed by the Navy Fleet 
Numerical Weather Central described in Haltiner (1971) but the FNWC 
analysis made only slight alterations to the field of x computed by the 
procedure above. It was felt that these alterations were not of suf- 
ficient importance to justify the additional program complexity and 
computer running time. 

If data were missing the station was rejected from consider- 
ation. This scheme, with slight input-output modifications, was found 
readily adaptable to derivation of temperature and dewpoint fields at 
all levels, although its ability at upper levels is hampered because data 


density is substantially reduced, 


— 
Acthiats. af 


D4) | 


_ ae ; © Vuh 
' % & 


: ah 7 ¥ 
r¢ 4, - —?7 —— = + 
1 Ty : at 2 e ‘ J a * 


+ bo 7 : 


aa i e! ® 
Ssert Es t. ieVioebks : r| : 3 i Sg. “ .i ae 
7 ® gid Li } be j f > 4 —" 7 
7 : 7 Ps 
. @ befvgmo ; ‘ omliapidin wettest 
; € a : P i 
— 
Fe aa 
“rug to 2en =) 2 A 
a ‘ 9efe>2 . a) 
'? 


SS - eee 7 
1 . 
_ 


~mhpasay apet Eatosi As cri Mela Ovi. 3a 


oe ee ; Ce 
Dine? saw yamine Fitts rugeua-soy~d 2) fala Oil 


CHAPTER 4 


OROGRAPHIC PRECIPITATION COMPUTATION 


4,1 Introduction 

The vertical velocities as calculated in this model may be a 
crude approximation of the real state of the atmosphere or at times even 
totally erroneous, Without doubt many important processes have been 
overlooked, For example, a study by Hess and Wagner (1948) of isentropes 
over the mountainous region of the northwestern United States clearly 
indicates that motion is much more complex than that depicted here, How- 
ever, there is at least one way of assessing the ascending givens 
computed, and that is through the realization of precipitation, It is 
not sufficient to evaluate simply the occurrence or non-occurrence of 
precipitation, but the recorded amounts of precipitation should bear some 
resemblance to the vertical motions responsible for the rain's formation. 

Attempts to calculate precipitation quantitatively have been 
made by Smagorinsky and Collins (1955) and by O'Neill (1966) based on 
computed vertical velocity fields. Both use a very simplified rain 
model which dates back to a paper by Fulks (1935). A similar approach 


is used here, 


4,2 Rate of Precipitation 


One of the most vital areas of meteorology and at the same time 


one of the most elusive is cloud physics, Ascending vertical motion is 
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a necessary but not sufficient condition for the formation of rain, 
Important refinements to the theory of the growth of a single isolated 
rain drop have been made by Fukuta and Walter (1970), but as yet a totally 
satisfactory explanation and model of the growth of rain droplets in a 
cloud is lacking. In view of this and other difficulties a simple 
approach is taken, 
To utilize Fulks expression for the rate of precipitation a 
number of drastic assumptions are required: 
(a) condensation nuclei are sufficiently numerous 
so that condensation commences the moment the 
saturation temperature is achieved, 
(b) on the scale of motion considered, temperature 
changes for an air parcel occur PRESETS EcRb ky 
(c) water storage within a cloud is ignored, 
(d) no evaporation of precipitation takes place 


between the cloud and the ground, 


The above is a summary of a discussion found in Pettersen (1956b, 
Vol. 2 ). In a more concise fashion, it is assumed that any condensation 
is the result of pseudo-adiabatic cooling, and condensation which occurs 
at the saturation temperature results in precipitation at the ground, 
Consider a volume of air of unit cross-section and height AZ, 
The total water vapour ity the volume is oqAZ where 9 is the density of 
dry air and q is the mixing ratio. The rate of precipitation must equal 
the rate of condensation from the assumptions stated, Let 
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be the rate of precipitation and note that dg/dt is the condensation 
which results from removal of moisture from the volume (and hence the 
negative sign). The vertical velocity may be expressed as w = dZ/at and 


if it is possible to neglect changes in the thickness AZ then: 


d 
Pie =. = 0a AZ m2 
i = (4-2) 
but 
q = 0,622e 
p-e 
and 
sie Oe AU ey (p-e)22 = ¢ d(p-e) (4-3) 
az (p=e)¢ dZ aZ 
where e = saturation vapour pressure. The vapour pressure e is 


usually very small in comparison to the total atmospheric pressure at 


any level so that: 


Ope _ Vas Vorsecee 

aZ oF Bae a 
where Rg = gas constant for dry air, 
and T = the absolute temperature, 
and 0 = density of the air, 


Equation (4-3) may be simplified by expanding the terms in the square 
brackets and sudstituting for dp/dZ from (4-4), 


dq 0.622 de Be ee ) 
a, es R,T \d2 RAT 


This in turn may be substituted into equation (4-2). Then: 


Labo ee OnGae ( 
cf Rat 


de, &2.5 
a7 ba R,T ee (4-5) 
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At saturation, all the water vapour terms can be calculated so that a 
knowledge of w yields the precipitation rate. 


The expression de/dZ can be evaluated by the use of the chain 


rule, i.€. 
de . de dt (4-6) 


It is assumed that condensation occurs at temperatures below the freezing 
point of water, and thus for all temperatures in the normal meteoro- 


logical range: 
e = 6,11 exp ee - 1) (4-7) 
ee N2 As 


from the Clausius-Clapeyron equation and: 


a (& = myLy e 
aT R*¥TC 
where R* = universal gas constant 
Ly = latent heat of evaporation of water at 273° Kelvin. 


To determine dT/dZ, Fulks used an equation derived by Brunt (1933) for 


the pseudo-adiabatic lapse rate, but a simpler expression is given by 


Hess (1959): 


i sel =) 
at seep (4-8) 
Osoe2ev. 
eld emacs) 


where c.. = heat capacity at constant pressure for dry air, 


Pp 


This ignores the heat capacity of the liquid water and the change of the 


latent heat of evaporation with temperature, If the atmosphere is div- 


a) 
—_ cohen Rint 
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By 


ided into a number of layers, then the total precipitation at the ground 


from a column of unit cross-section is: 


MS 


Ps (4-9) 


If 


vies at 

There are two ways of viewing the situation, If equation 
(4-5) is considered from the point of view of a parcel of thickness AZ, 
the rate of precipitation due to that volume will decrease as the temp- 
erature decreases, even though w may remain constant. Eventually, at 
some very low temperature, virtually all the moisture will have cond- 
ensed and the precipitation rate will be reduced to zero. However, if 
the parcel just considered is replaced by another saturated parcel with 
initial conditions equal to the previous one, then over a particular 
position the rate of precipitation for layer i remains the rey The 
second case is far more instructive, and multiplication by a time period 
A t+ results in the precipitation which should be realized at that pos- 
ition over the time interval, 

O'Neill assumed that saturated conditions were present when- 
ever ascending motions were calculated, in order to simplify his precip- 
itation calculations, but this is not necessary. Let the total ascent 
over a time interval At be: 

he OE TGs (4-9) 
If the parcel is initially unsaturated then part of the ascent (Z,) is 
required to reach saturation. Thus only At‘ of the total time incre- 


ment At is really amenable to the treatment of Fulks, An effective 


vertical velocity w' may be determined by: 
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and Z 


if 


eile w'At (4-10) 


The lift required for saturation is: 


Z = Toes ido (4-11) 
y Caray Ce. 
where T, = initial temperature of parcel 
Tg@o = initial dewpoint temperature of parcel 
Yq = lapse rate of temperature 
Ydq = lapse rate of dewpoint temperature 


It can be shown that to a close approximation (see e.g, Pettersen (1956b 


Vol. e2) 3 
Ya - ‘aq = 0.83°K per 100 meters. 
Therefore 
TOGnaAl 
Ee oes go (4-12) 


4.3 The Rain Model 

Equation (4-12) was used to calculate an effective vertical 
velocity over the time period At. The parameters of temperature, dew- 
point temperature, vertical velocity and heights of surfaces for which 
the previous fields were known for the entire grid were entered, 

The normal case of the 850-mb surface situated above ground 
level was dealt with in a straight forward manner, The layer between 
the surface and 850 mb was divided into two sub-layers, The lower half 
was assigned the vertical velocity computed for the surface, The temp- 
erature and dewpoint at the middle of the first sub-layer were computed 


from the gradients from the surface to 850 mb,. The next layer was the 
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upper half of the one just considered plus the lower half of the 850 

to 700-mb layer. The last layer includes the upper half of the 850 to 
700-mb layer the lower half of that from 700 to 500 mb, The attributes 
of the 850 and 700 mb levels were assigned to the last two layers, 
respectively. 

If the vertical velocity in a layer was negative or zero, the 
precipitation for that layer was set to zero. Next, the layers with 
positive vertical velocities of greater than 10 cm/sec were rejected and 
the precipitation set to an arbitrarily high value. Vertical motion in 
excess of 10 cm/sec was thought to exceed the assumption that precipi-~ 
tation is immediately realized at the ground, Sawyer (1956) and others 
recognize that such is not the case, and that a considerable time lag 
exists between the initial condensation and growth to rain oD size, 
Heavy precipitation is likely to be associated with these sizable ascent 
rates, but a more ‘Sophisticated approach is required to evaluate it. 

If 0 <w <10 cm/sec, then the ascent required for saturation 
was determined from equation (4-11), The saturation temperature was 
calculated at this level, and the effective ascent rate w’ for a 6 hour 
period determined, If w' was less than zero, saturation was not achieved 
and the layer’s precipitation set to zero. If all the above conditions 
were avoided, the precipitation rate was calculated using equation (4-5) 
and multiplied by At equal to 6 hours, By summing the contribution 
of the layers, the rain for a given grid position was determined, The 
amount first calculated in centimeters was converted to inches for con- 
parison with surface reports of precipitation, 

When the surface was above the 850-mb level, the computed 


surface vertical velocity was retained, but the temperature and dewpoint 
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were replaced by a linear interpolation between 850 and 700 mb, The 
calculations then proceeded as outlined above, 

The vertical velocity at 500 mb was set equal to zero as the 
orographic vertical velocity at that level was zero everywhere, except 
in a test case where a positive value of 1 cm/sec was observed. 

Obviously this model requires considerable refinement but as 
a first approximation it performed reasonably well. The results are 


presented in Chapter 5, 


CHAPTER 5 
RESULTS 


5S<. eintroduction 

Three synoptic situations are presented in this chapter. In 
each case an analyzed stream-function map is given for the 850-mb, 
700-mb and 500-mb levels. The computed vertical velocities are analyzed 
on separate charts at the 850-mb and 700-mb levels because the patterns 
derived at one-eighth grid length with the original terrain are much 
more complex than those found by Charette (1971). Orographic vertical 
velocities at 500 mb are in all cases zero but the 500-mb stream- 
function analysis is shown to illustrate the performance of the p cal- 
culation when data are scarce, In addition, the orographic vertical 
velocities associated with surface winds are illustrated in each case, 
Isotachs of vertical velocity are drawn at 5, 10 and 40 cm/sec while y 
is analyzed at 30-gpm increments. 

Included with the first example is a subjective analysis of 
the 850-mb constant pressure map for comparison with the y chart 
computed for the same level. 

The rain rap was applied to the second case only. In this 
instance, maps of computed orographic rain, actual reported rain, and 


the objective temperature analysis at the surface are shown in addition 


to those mentioned above. 
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5.2 October 14, 1960 - oozl 


5e2el The Synoptic Situation 

A deep low is centred 500 kilometers south of Kodiak Island, 
Alaska, while a weak ridge of high pressure dominates the circulation 
of Western Canada, In the course of the next 48 hours, the low tracks 
slowly northeast to the northern extremity of the Gulf of Alaska, pro- 
viding a sustained westerly flow across the British Columbia coast, 
The associated frontal wave is observed to occlude out of the low and 
move eastward across British Columbia, Lee cyclogenesis over north- 


eastern British Columbia is evident by 06Z on October 15. 


5.262 Stream-Function Analyses 

As mentioned earlier a subjective analysis (see pies (5-1).) 
is given of the 850-mb pressure surface in this case only, It should 
be mentioned that it was analyzed with considerable care to maintain 
contour gradients consistent with wind reports as determined by a geo- 
strophic wind scale. The original data were expressed in geopotential 
feet and the conversion is made from these units to gpm. Despite these 
precautions, the stream-function analysis (Figure (5-2)) was substan- 
tially different over the northeastern portion of the analyzed map, 
By drawing an intermediate contour, a closed high circulation is noted 
over the central Yukon which is not suggested by the subjective height 
analysis. In addition, the y analysis indicates a sharp trough in 


northwestern Alberta of which there is no indication on the subjective 
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analysis, The strong gradient along the left edge of the analysis 
is thought to be due to the lack of reports in that area, 

Of all the objective } analyses, the one of the 700-mb wind- 
field (Figure (5-3)) is probably the poorest. In this case, it is felt 
that the iterative procedure has been pushed beyond its capability. 

The ‘i gradient over northern British Columbia appears excessive and is 
believed due to the lack of reports in that area. As well, it can be 
seen that Annette and Fort Nelson are too widely separated for either 
to influence the other in the iterative procedure, and similarly 
Whitehorse and Prince George. Where the flow is uniform this may not 
be too critical, but the presence of both a trough and a ridge in close 
proximity in this data sparse region would necessitate more information 
for better resolution. Over the southern region of the map where data 
are more abundant no peculiarities appear, In fact it is of interest 
to note a diffluence of the streamlines in the vicinity of the Alberta 
Range of the Rocky Mountains. Also of interest is the presence of a 
trough over northeastern British Columbia, substantiating the feature 
atro50 MDs 

At the 500-mb level, the radius of influence was expanded to 
700 km. The results shown in Figure (5-4) are encouraging, The main 
problem area in the top left corner of the map is attributed to same 
problem as above, Again, troughing is present over north-central 
British Columbia. On the charts from which the wind data were extracted, 
the analyst failed to notice this feature at any of the levels analyzed, 
The trough provides the upper-air support for a minor wave on the arctic 


front which travelled rapidly eastward across northern Saskatchewan and 


Manitoba in the 12 hours following these maps, 
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On the next series of maps at 12Z it will be seen that the 
trough has also propagated rapidly eastward and is no longer within the 
analyzed area. The recognition of such minor features by the strean- 
function analysis could prove to be operationally beneficial. In this 
case, it proved to be advantageous over a height-contour analysis at 


all three levels, 


5.2.3 Vertical Velocity Fields 

The terrain=-induced vertical velocities at the surface 
(Figure (5-5)) illustrate a complex pattern of ascending and descending 
motions; however, some generalizations can be made, First, a large region 
of descent is observed along the eastern slopes of the Rocky Mountains 
south of 60°N, A parallel but narrower region of ascent lies along 
western slope. Vertical motions on the west coast reflect the irregular 
topography of that region. Associated with an offshore flow over the 
Alaska Panhandle are marked regions of descent, while to the south 
unorganized ascent appears over Vancouver Island, the western slopes of 
the Coastal Mountains and northwestern Oregon, 

At the 850-mb level (Figure (5-6)) the flow pattern is better 
organized, showing a region of pronounced descent to the lee of Vancouver 
Island with strong ascent along its western edge, and on the windward 
slope of the Coastal Mountains, Also, the Alaska Panhandle is under the 
influence of stronger ascending motion in contrast to the surface, 

At the 700-mb surface rapidly alternating areas of descent and 
ascent are evident (Figure (5-6)) over high terrain, The abnormally tight 
gradient over the Stikine Mountains mentioned previously makes this region 


one of the most active. Large areas of organized motion are not apparent, 
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FIGURE (5-1) 
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OCTOBER Tah1ge 


Subjective analysis of contour heights of the 850-mb 


pressure surface on October 14, 1960 - 00Z. Contours at 30-gpm inter- 


vals. 
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FIGURE (5-2) Modified stream-function analysis at 850-mb level on 
October 14, 1960 ~ 00Z, Contours drawn at 30-gpm intervals, 
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FIGURE (5-3) Modified stream-function analysis at the 700-mb level 
on October 14, 1960 - 00Z. Contours drawn at 30 gpm intervals, 
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FIGURE (5-4) Modified stream-function analysis at the 500-mb level 
on October 14, 1960 - 00Z. Contours drawn at 30-gpm intervals, 
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FIGURE (5-5) Surface vertical velocity induced by flow over sloping 
terrain on October 14, 1960 ~ 00Z. The heavy solid line separates 
regions of ascending and descending motion. Isotachs are drawn at 5, 
10 and 40 cm/sec. The letters A and D represent ascent and descent 


respectively. 
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FIGURE (5-6) Terrain-induced vertical velocities at the 850-mb level 

on October 14, 1960 - 00Z. Ascending and descending motions are marked 
by the letters A and D respectively and are separated by a heavy solid 
Jine, Isotachs at 5, 10 and 40 em/SeCe Shaded areas indicate surface 


below terrain. 
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FIGURE (5-7) 


on October 14, 1960 - 00Z. Ascent and descent are denoted by the 


letters A and D respectively and separated by heavy solid line, Iso- 
tachs at 5, 10 and 40 cm/sec, 


Terrain-induced vertical velocities at the 700=mb level 
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5.3 October 14,1960 = 127 


5e30el Synoptic situation 

This series of maps is just 12 hours later than those studied 
in the previous section. The low now occupies a position near the centre 
of the Gulf of Alaska and the associated frontal wave is analyzed to the 
northwest of the Queen Charlotte Islands, The circulation at the upper 


levels has increased markedly along the west coast. 


5.3.2 Stream-Function Analyses 

As a whole, the stream-function analyses in this series are 
handled much better than the earlier ones (see Figures (5-8), (5-9), 
(5-10)). The stronger winds ahead of the low are reflected by a tight- 
ening of the y gradient over the Alaska Panhandle, At the 850-mb a 
closed high is drawn just west of Jasper and a ridge over the same 
region is noted at the 700-mb level. This feature is not visible at the 
500-mb level and is therefore presumed to be a reflection of flow around 


the major barrier of the Alberta Range. 


523-3 Vertical Velocity Fields 

At the surface (Figure (5-11)), many of the features noted 12 
hours earlier are still present. The descending motion noted along the 
Alaska Panhandle is slightly more dominant as the offshore flow extends 
as far south as Vancouver Island. Slow descent covers much of Alberta 
except for marked descent near Pincher Creek (PC). 

At the 850-mb level, the main change has been the increase in 
vertical motions associated with strengthening winds over the Alaska 


Panhandle, Also of interest are the vertical motions near the Swan Hills 
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ob, 
in central Alberta, and over the Cypress Hills in southwestern 
Saskatchewan (see Figure (5-12)). 

Again, the 700-mb vertical velocity chart (Figure (5-10)) is 
distinguisned by its lack of organized motion, but to a lesser degree 
than the 00Z case, Descent to the lee of Vancouver Island is well 
marked, as is ascent along western slope of the Coastal Mountains in 
southwestern British Columbia. Ascending motions of similar magnitude 
to those shown on 00Z map along the western slope of the Rocky Mountains, 
from west of Jasper northward, are still visible while descent prevails 


to the lee of the same region, 


5.3.4 Surface Temperature Objective Analysis 

An illustration (see Figure (5-14)) of the derived temperature 
field at the surface using the scalar objective analysis discussed in 
Chapter 3 is included here, A total of 78 stations are used as input so 
that reasonable detail is possible, The dewpoint temperature analysis 
is not shown, but the pattern is similar, Since the temperatures at the 
surface are normally reported in degrees Fahrenheit, the values given on 
the analyzed map are also in Fahrenheit. However, temperatures are con- 
verted to degrees Kelvin for input to the thermodynamic calculations of 
the rain model, 

At upper levels where fewer than 20 stations for the same 
area are available for eee data, the fields of temperature and dew= 


point temperature are much smoother than those at the surface, as would 


be expected. 


5.3.5 Orographic Precipitation 


The analysis of computed precipitation due to orographic ascent 
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is shown in Figure (5-15) and the verification at 18Z of October 14, 
1960 in Figure (5-16). A first inspection of the computed chart would 
indicate that the calculated rainfall amounts are too high. However, 
Some very interesting features are present, such as the “hole” in the 
precipitation pattern over south-central British Columbia. Large 
orographic precipitation amounts are calculated along the West Coast 

as anticipated, and also substantial amounts along the windward slope of 
the Alberta Range, Naturally evaluation is difficult along the coastal 
area aS rain associated with the frontal system is also present, All 
that can be said is that the orographic effect appears to enhance the 
rainfall which would occur if the terrain were not present. Farther 
north along the Alaska Panhandle observed amounts are considerably less 
than those near the Queen Charlotte Islands and are of the en order 
of magnitude as the computed values, Also Sitka (SIT) reports no prec- 
ipitation for the previous 6 hours and none is calculated. 

The last few statements may seem to indicate that the precip- 
itation over the Alaska Panhandle is largely orographic in origin, but 
too much should not be read into these comparisons as much more data 
are needed to substantiate these findings. Another successful predi- 
ction is the precipitation at Fort Simpson (FS) in the Northwest Terr- 
itories. No frontal activity can be identified as the likely cause so 
that the precipitation may be entirely orographic in nature. 

Computed trace amounts should be regarded with suspicion as 
the assumption that no evaporation occurs during the descent is not 
strictly valid. Another problem which is not treated at all by the rain 
model is the horizontal advection of raindrops. O'Neill (1966) found 


that computed maxima were displaced upstream from the observed precip- 


ot sain) to of « : honeatse ieee) i+, ry (etney : , 


kines gaa losugtios erty te igoitowren? L2Sk 


_~feveaccX ist? of Gn gfavtie Lleinivr Gegee 


Yt en x Foti £ () Ac cyedie iis! rit @ ores 
; eed a » = =5 ee i 
‘ ; : ; ; 
fracdr-s an oc g GIOAIuOlas ot =p — sey t Dige aD 
«4 
i a, jy Thi nl spleens “ipsa toe oO "7 
‘ +f j ‘p44 a - _< 


. i ,> a 
Leg eh 483 ah) dohisotees AL teal Per ated 


‘ 7 " , A es on Ze ides as - 
‘ has, 4 j — ri ae $£ we 7 —_ i . ai 
z yy = 7 > 
“r pers “7 : in 2} “7 rent w mT un: Lad Aa we: pea 5 8 nao FP tehee 


= . 


fips ay tad Tl aesofilvew get Lie 
iota ay ns niithel giaeth ads piots 
at ad seh septs 7 iz 
atau On pent od stvVeT? C62 Mi B 

_ te asael i ited ‘e- dah cael si? 
| wind batter : : wit lah we cave okies ¢ 
- aa x 
: : Th sid aie «6 ego! fr sgaredue Se on “asl 


~—qeT fepetirraiwat? tot ) teat =—=t? Ja epomeares cits 


=% 1% bet 


on Ge Viet. Jac cima > 


® 3n 2a 2 


dakia om lens 
> ay 


+S 


: “ 
rea aor Ge ae 
tall 


2 


itation maxima, but no statement concerning this matter can be made 
here as observations are far too scarce, 

One of the worst failures of the model was a computed value of 
.05 inches at Pincher Creek (PC), Alberta while the actual 18Z report 
indicates virtually clear skies and no precipitation over the past 6 
hours. The explanation is believed to lie in the assumption that the 
circulation did not change markedly over the 6-hour period. 

The lack of surface observations makes further comparisons 
difficult. A better evaluation could be attempted by supplementing 
regular reports with climatological reports and forestry observations in 
the summer months. However, the above discussion does offer some encour- 


agement that the results may be useful. 
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ied stream-function analysis at the 850-mb level on 


FIGURE (5-8) Modif 


Contours are drawn at 30-gpm intervals. 


October 14, 1960 - 122. 


_ _ 7 - — 
° » < 
4 7 
_—_"* _ - 
. =, 
i) 
~*~, 
{ —_ 
4 a 
' i. 
8 
/ 
\ | 
i 
\ 
4 
\ 
* 
i 
} 
| 
| 
a 
J ; 
j 
, 7 
j 
ce | 
* 
A 
_ 
4, 
‘ . q P 
| ' ol 
: . ¥ ss 
i. ae 
. . - =y ad * 


p Dagens wher, a¢ 
7 & ‘ 


ee) | reese 


a 7 7 7 _ : 7 a? 
Bese dp clafisee ssh ros/ pa 


Say. ae ean °c @& ro, 
' a4 


OPr~. : 
fi E. yh 
OPIOBER HH 1200 12Z 


/ 


FIGURE (5-9) Modified stream-function analysis at the 700-mb level 
on October 14, 1960 - 12Z. Contours drawn at 30-gpm intervals, 
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FIGURE (5-10) Modified stream-function analysis at the 500-mb level 
Contours drawn at 30-gpm intervals. 


on October 14, 1960 ~ 12Z, 


FIGURE (5-11) Surface vertical velocity induced by flow over sloping 
terrain on October 14, 1960 - 12Z. The heavy solid line separates 
regions of ascending and descending motion. Isotachs are drawn at 5, 
10 and 40 cm/sec. The letters A and D represent ascent and descent 
respectively. 
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FIGURE (5-12) Terrain-induced vertical velocities at the 850-mb level 
on October 14, 1960 - 12Z. Ascending and descending motions are marked 
by the letters A and D respectively and are separated by a heavy solid 
line, Isotachs at 5, 10 and 40 cm/sec. Shaded areas indicate surface 
below terrain, 
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FIGURE (5-13) Terrain-induced vertical velocities at the 700-mb level 


on October 14, 1960 - 12Z. Ascent and descent are denoted by the 
letters A and D respectively and separated by heavy solid line. 
tachs at 5, 10 and 40 cm/sec. 
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FIGURE (5-14) Objective analysis of surface temperature on October 14, 
1960 - 122. Isotherms are drawn every 5 degrees and temperatures are 
in Fahrenheit. 
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FIGURE (5-15) Analysis of computed rainfall amounts in inches, 
Isohyets are drawn at .00 (a trace), .05, .10, and .15 inches. 
denote ascent greater than 10 cm/sec above the point considered, 
are calculated for October 14, 1960 from 12Z to 182. 
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FIGURE (5-16) Observed 6-hour rainfall amounts over past 6 hours on 
October 14, 1960 - 18Z. Also shown is present weather, Each station 
used in the analysis is marked by a large dot. The presence of an 
amount to the lower right of a dot gives a value greater than zero for 
the precipitation during the previous 6 hours, 
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5.4 June 21, 1960 - 002 


54.1 The Synoptic Situation 

This particular case is chosen to illustrate the point that 
the Rocky Mountains are partially responsible for much of the precip- 
itation received in Alberta in the spring and early summer months, for, 
as Reinelt (1970) has pointed out, an easterly circulation for Alberta 
is just as much an upslope situation as a westerly flow is for the 
British Columbia coast, 

At 00Z, June 19, 1960, a somewhat ill-defined low and an 
associated frontal wave are situated over southern British Columbia, 
The associated low at 500 mb is centred near the southern tip of the 
Queen Charlotte Islands, By the same time on June 20, two surface 
centres are visible, one south of Edmonton and the other immediately 
west of Glasgow, Montana, The 500-mb low has now assumed the former 
position of the surface feature. The temperature at its centre is near 
~30°C, 

By June 21 the Glasgow low has merged with a series of lows 
which comprise a trough complex extending from Colorado to southeast 
Saskatchewan, recurving to the low some 100 kilometers north of Edmonton, 
The cold low is now established almost vertically above this surface 
position, and the air-mass associated with it is moist and unstable, 

The surface low eventually intensified and moved slowly south- 
east over the next two days, during which time much of Alberta was sub- 


jected to a northeasterly flow with its consequent cloud and showers, 


5.4.2 The Stream-Function Analyses 


All three levels are handled well by the objective stream- 
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function analysis (see Figures (5-17), (5-18), (5-19)). A peculiar 
gradient is evident in the upper right corner of the 850-mb analysis but 
the absence of a reported wind at Inuvik is a likely cause. Also, a 
direct iterative compromise between Norman Wells and Fairbanks is not 
possible as they are separated by a ae ence greater than 600 kn, 

The only feature not picked up by the subjective height- 
contour analysis which was recognized by the ‘ analysis was a minor 


trough at the 850-mb level as indicated by the wind at Calgary, 


5.4.3 Vertical Velocity Fields 

The instability of the air mass which covers the southern half 
of Alberta is the reason for much of the precipitation observed in this 
area. However, the most persistent region of rainfall is observed over 
the Peace River een and south to Jasper. With reference to vertical 
velocity charts, Figures (5-20) to (5-22), this same region is marked 
by substantial orographic ascent at all three levels. The position of 
the low makes it impossible for ascent to occur along the entire eastern 
slope of the Alberta Range. In fact, an area of descent is observed to 
extend almost unbroken from Edson south where a return westerly flow 
exists, 

In this instance, the rain is not solely an orographic phen- 
omenon but, doubtless, terrain-induced vertical motions play an import- 
eigen & SOC. to USL exactly what this role may be is masked by the 
presence of convective precipitation. The rain model is not capable of 


handling instability cases so that it was not applied to this situation. 
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FIGURE (5-17) Modified stream-function 
June 21, 1960 - 00Z. Contours drawn at 


analysis at 850-mb level on 
30-gpm intervals, 
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FIGURE (5-18) Modified stream-function 
on June 21, 1960 - 00Z. Contours drawn 
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analysis at the 700=-mb level 
at 30-gpm intervals. 
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FIGURE (5-19) 


Modified stream-function analysis at the 500=-mb level 


on June 21, 1960 - 00Z,. Contours drawn at 30-gpm intervals, 
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FIGURE (5-20) Surface vertical velocity induced by flow over sloping 
terrain on June 21, 1960 - 00Z. The heavy solid line separates 
regions of ascending and descending motion. Isotachs are drawn at 5, 
10 and 40 cm/sec. The letters A and D represent ascent and descent 
respectively. 
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FIGURE (5-21) Terrain-induced vertical velocities at the 850-mb level 
on June 21, 1960 - 00Z. Ascending and descending motions are marked 
by the letters A and D respectively and are separated by a heavy solid 


line, Isotachs at 5, 10 and 40 cm/sec, Shaded areas indicate that 
this level is below the terrain, 
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FIGURE (5-22) Terrain-induced vertical velocities at the 700-mb level 
on June 21, 1960 - 00Z. Ascent and descent are denoted by the letters 
A and D respectively and separated by heavy solid line. Isotachs at 
5, 10 and 40 cm/sec. 
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CHAPTER 6 
CONCLUSIONS 


The iterative procedure for calculating the modified strean- 
function has been improved to the point that far fewer station reports 
are required to arrive at a reasonable field, After considerable exper- 
imentation with Charette's (1971) model, this approach was found neces- 
sary because a great deal of time was spent subjectively interpolating 
Winds and heights at the particular pressure level under examination from 
charts already subjectively analyzed by hand, This amounted to placing 
a great deal of faith in the geostrophic relation. Therefore, it was 
decided to revert to Endlich and Mancuso's original stream-function 
calculation which was based solely on the reported winds, The reported 
heights at radiosonde stations were used only as an initial guess for 
the modified stream-function in order that the final field had values 
Closely approximating the heights at the constant pressure level. 

This approach proved reasonably satisfactory. The regular 
radiosonde reports supplemented by pilot balloon reports were sufficient 
to determine the flow over the continent, Because of the lack of ade- 
quate data networks this is not possible over the oceans so that this 
procedure is restricted primarily to the land areas. Emphasis should 
be placed on the need for the supplementary reports, particularly in 
northern British Columbia and the central Yukon, As is illustrated 


in the 00Z, October 14, 1960 situation, some unusually strong gradients 
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may arise when the distance between stations is large and the inter- 
vening pattern not uniforn, 

The continuous weighting scheme introduced here, while only 
Slightly less arbitrary, substantially simplifies the computer pro- 
gramming and reduces the storage required, It also has the advantage 
that minor peculiarities may be maintained in those regions where data 
are dense, while providing a reasonably smooth interpolation in data- 
sparse areas, 

The final stream-function objective analysis provides a 
useful tool and could easily be adapted to a small computer facility. 
An error rejection scheme would have to be formulated for direct input 
of information from teletype. In this study, analysis was restricted 
to the standard levels of synoptic analysis, but Speretionelly it could 
be applied to any level where reasonable data density exists, 

Little change was made to the method of computing orographic 
vertical velocities proposed by Charette, The use of the extracted 
terrain heights at one-eighth grid length intervals combined with the 
reduction of the finite difference increment for calculation of terrain 
Slopes complicated the pattern of orographic vertical motion. It is 
difficult to assess whether this yields any significant improvement or 
merely confuses the issue, Because of friction, vertical velocities at 
the surface are likely. to be underestimates, while at the 850-mb and 
700-mb surfaces the estimates may be too high, Daily evaluation is 
required to reveal systematic errors so that a more realistic decrease 
of the orographic influence with height may be formulated. 

| A precipitation model similar to the one suggested in Chapter 4 


is useful in the evaluation of ascending motion but provices no prac- 
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tical estimate of the magnitude of downward motion, Charette's treatment 
of valleys and peaks may provide a better method for narrow valleys, 
but to the lee of major topographic features would tend to underestimate 
negative vertical velocity. The current procedure is probably better 
in the latter case, provided one realizes that minor centres of descent 
over frequently fluctuating topography are not likely present, due to 
stagnation of air in valleys. 

As suggested in Chapter 5, a finer network of rain gauges is 
required to evaluate the orographic precipitation. Climatological and 


forestry reports would be invaluable as well, 
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APPENDIX A 


Calculation of Wind Direction With Respect to the Grid 


(180°- ) 
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BO -(18C0 =A) + (18020) = 10° 


FIGURE (A-1) Angle the locally measured wind makes with the x-axis of 
the grid, 


With reference to Figure (A-1), it may be seen that the wind 
direction with respect to the grid depends on the local wind direction 
D, which is measured relative to true north, the longitude A where the 
wind is observed, and the orientation of the grid. The grid is rotated 
10° so that the ordinate is parallel to 80°W and the abscissa to 


10°R, Therefore, once the wind direction is determined with respect to 


ae) 


the zero meridian, subtraction of ten degrees yields the wind direction 


with respect to the abscissa of the grid. 


The angle the wind makes with the local meridian is 180° - D 


while the supplement of the longitude is 1809 - , 


The direction of the wind relative to the x-axis of the grid 


ASB 
(e0CF=9) + (180° = 91) F-71083 5098— D — 9h, (A-1) 


from which the wind components may be calculated, 
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